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Summary. We consider random Schrodinger equations on or 7 /^ for d > 3 with 
uncorrelated, identically distributed random potential. Denote by A the coupling 
constant and i/)f the solution with initial data ?/)o. Suppose that the space and time 
variables scale as a; ~ t ~ with 0 < k < /to, where kq is a sufficiently 

small universal constant. We prove that the expectation value of the Wigner distri¬ 
bution of ipt, EVU/,j (a;, n), converges weakly to a solution of a heat equation in the 
space variable x for arbitrary initial data in the weak coupling limit A —> 0. The 
diffusion coefficient is uniquely determined by the kinetic energy associated to the 
momentum v. 


1 Introduction 

Brown observed almost two centuries ago that the motion of a pollen suspended 
in water was erratic. This led to the kinetic explanation by Einstein in 1905 that 
Brownian motion was created by the constant “kicks” on the relatively heavy pollen 
by the light water molecules. Einstein’s theory, based upon Newtonian dynamics 
of the particles, in fact postulated the emergence of the Brownian motion from a 
classical non-dissipative reversible dynamics. Einstein’s theory became universally 
accepted after the experimental verification by Perrin in 1908, but it was far from 
being mathematically rigorous. 

The key difficulty is similar to the justihcation of Boltzmann’s molecular chaos 
assumption (Stofizahlansatz) standing behind Boltzmann’s derivation of the Boltz¬ 
mann equation. The point is that the dissipative character emerges only in a scaling 
limit, as the number of degrees of freedom goes to infinity. 

The hrst mathematical definition of the Brownian motion was given in 1923 by 
Wiener, who constructed the Brownian motion as a scaling limit of random walks. 
This construction was built upon a stochastic microscopic dynamics which by itself 
are dissipative. 

The derivation of the Brownian motion from a Hamiltonian dynamics was not 
seriously investigated until the end of the seventies, when several results came out 



2 


Laszlo Erdos, Manfred Salmhofer, and Horng-Tzer Yau 


almost simultaneously. Kesten and Papanicolaou [16] proved that the velocity dis¬ 
tribution of a particle moving in a random scatterer environment (so-called Lorenz 
gas with random scatterers) converges to the Brownian motion in a weak coupling 
limit for d > 3. The same result was obtained in d = 2 dimensions by Diirr, Gold¬ 
stein and Lebowitz [8]. In this model the bath of light particles is replaced with 
random static impurities. In a very recent work [18], Komorowski and Ryzhik have 
controlled the same evolution on a longer time scale and proved the convergence to 
Brownian motion of the position process as well. 

Bunimovich-Sinai [5] proved the convergence of the periodic Lorenz gas with a 
hard core interaction to a Brownian motion. In this model the only source of ran¬ 
domness is the distribution of the initial condition. Finally, Diirr-Goldstein-Lebowitz 
[7] proved that the velocity process of a heavy particle in a light ideal gas converges 
to the Ornstein-Uhlenbeck process that is a version of the Brownian motion. This 
model is the closest to the one in Einstein’s kinetic argument. 

An analogous development happened around the same time towards the rigorous 
derivation of the Boltzmann equation. It was proved by Gallavotti [14], Spohn [26] 
and Boldrighini, Bunimovich and Sinai [3] that the dynamics of the Lorenz gas with 
random scatterers converges to the linear Boltzmann equation at low density on the 
kinetic time scale. Lanford [19] has proved that a truly many-body classical system, 
a low density gas with hard-core interaction, converges to the nonlinear Boltzmann 
equation for short macroscopic times. 

Brownian motion was discovered and theorized in the context of classical dynam¬ 
ics. Since it postulates a microscopic Newtonian model for atoms and molecules, it 
is natural to replace the Newtonian dynamics with the Schrodinger dynamics and 
investigate if Brownian motion correctly describes the motion of a quantum particle 
in a random environment as well. One may of course take first the semiclassical limit, 
reduce the problem to the classical dynamics and then consider the scaling limit. 
This argument, however, does not apply to particles (or Lorenz scatterers) of size 
comparable with the Planck scale. It is physically more realistic and technically con¬ 
siderably more challenging to investigate the scaling limit of the quantum dynamics 
directly without any semiclassical limit. We shall prove that Brownian motion also 
describes the motion of a quantum particle in this situation. It is remarkable that 
the Schrodinger evolution, which is time reversible and describes wave phenomena, 
converges to a Brownian motion. 

The random Schrodinger equation, or the quantum Lorentz model, is given by 
the evolution equation: 

idtMx)^HMx), H ^ + Wu^ix) ( 1 ) 

where A > 0 is the coupling constant and VL is the random potential. 

The first time scale with a non-trivial limiting dynamics is the weak coupling 
limit, A ^ 0, where space and time are subject to kinetic scaling and the coupling 
constant scales as 

t ^ te~^, X ^ xe~^, A =-x/i . (2) 

Under this limit, the appropriately rescaled phase space density (Wigner distribu¬ 
tion, see (10) later) of the solution to the Schrodinger evolution (1) converges weakly 
to a linear Boltzmann equation. This was first established by Spohn (1977) [25] if 
the random potential is a Gaussian random field and the macroscopic time is small. 
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This method was extended to study higher order correlations by Ho, Landau and 
Wilkins [13]. A different method was developed in [10] where the short time restric¬ 
tion was removed. This method was also extended to the phonon case in [9] and to 
the lattice case in [6]. 

For longer time scales, one expects a diffusive dynamics since the long time limit 
of a Boltzmann equation is a heat equation. We shall therefore take a time scale 
longer than in the weak coupling limit (2), i.e. we set t ~ n > 0. Our aim 

is to prove that the limiting dynamics of the Schrodinger evolution in a random 
potential under this scaling is governed by a heat equation. This problem requires 
to control the Schrodinger dynamics up to a time scale This is a much harder 

task than first deriving the Boltzmann equation from Schrodinger dynamics on the 
kinetic scale and then showing that Boltzmann equation converges to a diffusive 
equation under a different limiting procedure. Quantum correlations that are small 
on the kinetic scale and are neglected in the first limit, may contribute on the longer 
time scale. 

We consider two models in parallel. In the discrete setup we put the Schrodinger 
equation (1) on Z'*, i.e. we work with the Anderson model [2]. Thus the kinetic 
energy operator on £^(Z'^) is given by 

{Af){x) :=2d f{x)- ^ f{x + e) (3) 

|e| = l 


and the random potential is given by 

Vij{x) = Vj{x) , V~f{x) := v~fS{x — ■y) (4) 

7ez<i 

where v^y are real i.i.d. random variables and 5 is the lattice delta function, 5(0) = 1 
and 5{y) = 0, y / 0. 

In the continuum model we consider the usual Laplacian, as the kinetic 

energy operator on L^(R‘*). The random potential is given by 

VL(a;)= [ B{x-y)dfi^{y), (5) 

where is a Poisson process {y-y : 7 = 1, 2,. ..} on R'^ with unit density and i.i.d. 
random masses, Vj, i.e. Vj5{- — y-y), and B : R'* ^ R is a smooth, radially 

symmetric function with rapid decay, with 0 in the support of B. 

Since we investigate large distance phenomena, there should be no physical dif¬ 
ference between the continuum and discrete models. On the technical level, the 
discrete model is more complicated due to the non-convexity of the energy surfaces 
of the discrete Laplacian in momentum space. However, the continuum model also 
has an additional technical difficulty: the large momentum regime needs a separate 
treatment. 

Our proof builds upon the method initiated in [10]. In that paper the continuum 
model with a Gaussian random field was considered. Here we also consider the 
discrete model and non-Gaussian randomness, in order to demonstrate that these 
restrictions are not essential. On the Boltzmann scale this extension has also been 
achieved by Chen [6]. The other reason for working on the lattice as well is to make 
a connection with the extended state conjecture in the Anderson model. 
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We recall that the Anderson model was invented to describe the electric conduc¬ 
tion properties of disordered metals. It was postulated by Anderson that for localized 
initial data the wave functions for large time are localized for large coupling constant 
A and are extended for small coupling constant (away from the band edges and in 
dimension d > 3). The localization conjecture was hrst established rigorously by 
Goldsheid, Molchanov and Pastur [15] in one dimension, by Frohlich-Spencer [12], 
and later by Aizenman-Molchanov [1] in several dimensions, and many other works 
have since contributed to this held. The extended state conjecture, however, has 
remained a difficult open problem and only very limited progress has been made. 

Most approaches on extended states focused on the spectral property of the 
random Hamiltonian. It was proved by Klein [17] that all eigenfuctions are extended 
on the Bethe lattice. In Euclidean space, Schlag, Shubin and Wolff [24] proved that 
the eigenfunctions cannot be localized in a region smaller than for some 5 > 0 

in d = 2. Chen [ 6 ], extending the method of [10] to the lattice case, proved that the 
eigenfunctions cannot be localized in a region smaller than A“^ in any dimension 
d > 2 with logarithmic corrections. Lukkarinen and Spohn [21] have employed a 
similar technique for studying energy transport in a harmonic crystal with weakly 
perturbed random masses. 

A special class of random Schrodinger equation was proposed to understand the 
dynamics in the extended region. Instead of random potential with i.i.d. random 
variables, one considers a random potential VL(®) with a power law decay, i.e., 

Vui{x) = h{x)uJx , h{x) ~ jxp’’ 

where are mean zero i.i.d. random variables and > 0 is a fixed parameter. 

If ?; > la standard scattering argument yields that for A small enough 
has absolutely continuous spectrum. Using cancellation properties of the random 
potential, Rodnianski and Schlag [22] have improved the same result to 77 > 3/4 in 
d > 2 and recently, J. Bourgain [4] has extended it to 77 > 1/2. Eor 77 > 1/2 the 
particle becomes essentially ballistic at large distances and there are only finitely 
many effective collisions. 

In summary, in all known results [24, 22, 4, 6 ] for the Anderson model (or its 
modihcation) in Euclidean space the number of effective collisions are hnite. In the 
scaling of the current work (13), the number of effective scatterings goes to infinity 
in the scaling limit, as it should be the case if we aim to obtain a Brownian motion. 

As in [ 6 ], our dynamical result also implies that the eigenfunctions cannot be 
localized in a region smaller than for some <5 > 0 and dimension d > 3 (one 

can choose 5 = k/ 2 with k from Theorem 1). Though this result is the strongest in 
the direction of eigenfunction delocalization, we do not focus on it here. 

Our main result is that the time reversible Schrodinger evolution with random 
impurities on a time scale is described by a dissipative dynamics. In fact, this 

work is the hrst rigorous result where a heat equation is established from a time 
dependent quantum dynamics without hrst passing through a semiclassical limit. 

In this contribution we explain the result and the key ideas in an informal man¬ 
ner. The complete proof is given in [11]. 
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2 Statement of main result 

We consider the discrete and the continuum models in paralell, therefore we work 
either on the d-dimensional lattice, Z'*, or on the continuous space, R'*. We always 
assume d > 3. Let 

:=AK; (6) 

denote a random Schrodinger operator acting on 7d = or Ti. = L^(R‘*). The 

kinetic energy operator and the random potential are defined in (3)-(5). We assume 
that Ev-y = Eu® = 0, E«^ = 1 and Eu^‘* < oo. 

In the discrete case, the Fourier transform is given by 

Tip) = (^/)(p) := ^ , 

where p — {p^^\ ■ ■ ■ ,p^'^^) € := [—i, i]'*. Sometimes an integral notation will be 

used for the normalized summation over any lattice (dZ)'*: 

f (■■■)dx ■- d'* ^ (•••)• 

xe(sz)<i 

The inverse Fourier transform is given by 

In the continuous case the Fourier transform and its inverse are given by 

i^f )ip)~[ f{x)dx , {T~^g){x)=( . 

We will discuss the two cases in parallel, in particular we will use the unified integral 
notations J (• • •)d® and J (• • •)dp. The letters x, y, z will always be used for position 
space coordinates (hence elements of (dZ)‘* or R'*). The letters p, q, r, u, v, w denote 
for d-dimensional momentum variables (elements of (T/d)‘* or R'*). 

The Fourier transform of the kinetic energy operator is given by 

/) (P) = e(p)/(p) . 

The dispersion law, e(p), is given by 
d 

e(p) := ^~~^(1 — cos(27rp^'^)), and e{p) := —p^ 

i=l 

in the discrete and in the continuous case, respectively. 

For h : T'* ^ C and an energy value e £ [0, 2d] we introduce the notation 

[h]{e) ■-J h{v)S{e - e{v))dv ■.= J h{q) (7) 
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where du{q) — dv^iq) is the restriction of the d-dimensional Lebesgue measure to 
the level surface Ee := {q : e{q) = e} C T^*. By the co-area formula it holds that 

[h]{e)de = J h{v)dv . (8) 

We define the projection onto the energy space of the free Laplacian by 

He) ■■= Wie) = j S{e - e{u))du . (9) 

In the continuous case we define analogous formulas for any function h : R'* —^ C 
and energy value e > 0. 

Define the Wigner transform of a function ip £ or 'll; € L^(R‘*) via its 

Fourier transform by 

W^{x,v) := J . 

In the lattice case the integration domain is the double torus (2T)‘* and x runs over 
the refined lattice, x £ (7^/2)'^. For e > 0 define the rescaled Wigner distribution as 

W^X, V) := (^, f) . (10) 

(with X £ {e7/2)'^ in the lattice case). 

The weak coupling limit is defined by the following scaling: 

T := et, X := ex, e = )? . (11) 

In the limit e ^ 0 the Wigner distribution Wfj, (XtV) converges weakly to the 
Boltzmann equation ([10], [6]) 

(dr + ^ Ve(F) • Xx^ Fr{X, V) = j dUa{U, V) Ft{X, U) - Ft{X, V) (12) 

where ^Ve(F) is the velocity. The collision kernel is given by 
cr([/, V) := 2iT5{e{U) — e(F)) discrete case 

(j{U, V) := 27r|i3(f7 — V)\^6{e{U) — e{V)) continuous case . 

Note that the Boltzmann equation can be viewed as the generator of a Markovian 
semigroup on phase space. In particular, the validity of the Boltzmann equation 
shows that all correlation effects become negligible in this scaling limit. 

Now we consider the long time scaling, i.e. with some «: > 0, 



a; = t = A“'"“V = e = 


(13) 
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Theorem 1. [Quantum Diffusion on Lattice] Let d = 3 and ipo £ be an 

initial wave function with V'o € C'^(T'^). Let ipit) = solve the Schrodinger 

equation (1). Let 0{x,v) be a function on R'* x T'* whose Fourier transform in the 
first variable, denoted by 0{^,v), is a function on R'^ x T'* and 



du|(!?(C,u)||^| <C. 


Fix e e [0,2d]. Let f be the solution to the heat equation 

drfiT, X, e)=Xx- D{e)Vxf{T, X, e) 
with the initial eondition 


f{0,X,e) 


5(X) (e) 


(14) 


(15) 


and the diffusion matrix D 

( sin(27ru*'d) • sin(27ru^'’^) \ 

^ hi = 1,2,3. (16) 

Then for k < 1/2000 and e and \ related by (13), the Wigner distribution satisfies 


lim 

e^O 


(eZF 


dX / dvd{X,v)EW;^^-^-.j,^iX,v) 




(17) 



dvd{X,v)f{T,X,eiv)) . 


By the symmetry of the measure (•)e under each sign flip Vj —vj we see that 
D{e) is a constant times the identity matrix: 


Dij{e) = De Si. 


D, := 


^ sin^(27rw *'^^)') 
2tv F{e) 


in particular we see that the diffusion is nondegenerate. 

The diffusion matrix can also be obtained from the long time limit of the Boltz¬ 
mann equation (12). For any fixed energy e, let 


Lef{v) := / dtt a{u,v)[f{u) - f{v)], e{v) = e 


(18) 


be the generator of the momentum jump process on Xe with the uniform stationary 
measure {■}(,■ The diffusion matrix in general is given by the velocity autocorrelation 
function 

poo 

Dij{e) = / dt ^ sin(27rud)(t)j . sin(27ru*''’l(0)))^ , (19) 

Jo 

where v{t) is the process generated by Le. Since the collision kernel g{U,V) is 
uniform, the correlation between v(t) and u(0) vanishes after the first jump and we 
obtain (16), using 
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J dw a(u,v) = 2'K<P{e) , e{v) = e . 

The result in the continuum case is analogous. The diffusion matrix is again a 
constant times the identity matrix, Dij{e) — and Df, is again given by the 

velocity autocorrelation function 

De := (u(t) • u(0))e (20) 

using the spatial isotropy. In this case De cannot be computed as a simple integral 
since the outgoing velocity u in the transition kernel a{u,v) of the momentum 
process depends on the direction of v. 

Theorem 2. [Quantum Diffusion on Let d = 3 and ipo G L^(R‘^) be an initial 
wave function with |■^o('l’)P|w|^ € for a sufficiently large N. 

Let ip{t) = 'ipt,w solve the Schrodinger equation (1). Let 0{x,v) be a function 
whose Fourier transform in x, denoted by 0{£,,v), is a function on R'* x R'* and 


dCdu|0(C,v)||C| <C7. (21) 

Let e > 0 and let f be the solution to the heat equation 

drfiT, X, e) = De Axf{T, X, e) (22) 

with diffusion constant given in (20) and with the initial condition 

f{0,X,e) :=5(Y)[|V^o(«)|"](e). 

Then for k < 1/500 and e and A related by (13), the Wigner distribution satisfies 


lim 


R'ixR'^ 


AXAvO{X, v)EW;^^.^- 2 r) {X, v) 


(23) 


dXdv 0{X,v)f{T, X,e(v)) 


R'^xR'i 


The main tool of our proof is to use the Duhamel expansion to decompose the 
wave function into elementary wavefuctions characterized by their collision histories 
with the random obstacles. Assume for the moment that the randomness is Gaus¬ 
sian and high order expectations can be computed by Wick pairing. The higher 
order cumulants arising from a non-Gaussian randomness turn out to be negligible 
by a separate argument. Therefore, when computing the expectation of a product 
involving ip and ip (e.g. E W,p), we pair the obstacles in the collision histories of ip 
and Ip and we thus generate Feynman graphs. 

If we take only the Laplacian as the free part in the expansion, even the am¬ 
plitudes of individual graphs diverge in the limit we consider. However, this can be 
remedied by a simple resummation of all two-legged insertions caused by the lowest 
order self-energy contribution The resummation is performed by choosing an ap¬ 
propriate reference Hamiltonian Hq for the expansion. After this rearrangement, all 
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graphs have a finite amplitude in our scaling limit, and the so-called ladder graphs 
give the leading contribution. 

Each non-ladder graph has a vanishing amplitude as A ^ 0 due to oscillatory 
integrals, in contrast to the ladder graphs where no oscillation is present. However, 
the number of non-ladder graphs grows as k\, where k ~ ~ A”*^ is the typical 

number of collisions. To beat this factorial growth, we need to give a very sharp 
bound on the individual graphs. 

We give a classification of arbitrary large graphs, based on counting the number 
of vertices carrying oscillatory effects. The number of these vertices is called the 
degree of the graph. For the ladder graphs, the degree is zero. For general graphs, 
the degree is roughly the number of vertices after removing all ladder and anti-ladder 
subgraphs. We thus obtain an extra A'^ power (for some c > 0) per non-(anti)ladder 
vertex. This strong improvement is sufficient to beat the growth of the combinatorics 
in the time scale we consider. To our knowledge, nothing like this has been done in 
a graphical expansion before. 

For a comparison, the unperturbed Green functions in the perturbation ex¬ 
pansion for the many-fermion systems for small temperature and for the random 
Schrodinger equation for large time are given by 

1 1 
*Po + — A* ’ p"^ — a + ir] 

In the many-fermion case, po £ Mp = {^(2n -I- 1) : n £ Z} where (3 ~ T~^ is the 
inverse temperature. In the random Schrodinger case, rj ^ . Their properties 

are different: 

i ^ y dp\ipQ+p‘^ - p 

PoGMp 

Notice the divergence is more severe for the random Schrodinger equation case. 

Finally we note that the threshold k < 1/2000 in our theorem can be significantly 
improved with more detailed arguments. However, one cannot go beyond k — 2 with 
only improvements on estimates of the individual graphs. The Duhamel formula 
must be expanded at least up to fc = A^t = A“'‘, which is the typical number 
of collisions up to time t. Even if one proves for most graphs the best possible 
estimate, it cannot beat the k\ combinatorics when k ^ A“^, i.e., ^ 1 

for k ^ A“^. A different resummation procedure is needed beyond this threshold to 
exploit cancellations among these graphs. 


|log/3|, 


p — a + iri\ 


3 Sketch of the proof 

We present the main ideas of the proof for the lattice case and comment on the 
modifications for the continuous case. 

3.1 Renormalization 

Before expanding the solution of the Schrodinger equation (1) via the Duhamel 
formula, we perform a renormalization of the ” one-particle propagator” by splitting 
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the Hamiltonian as H — Hq + V, with Hq already containing the part of the self¬ 
energy produced by immediate recollisions with the same obstacle. This effectively 
resums all such immediate recollisions. 

Let 6{p) := 0{e{p)), where 0{a) := lime^o+ 0s(a) and 

0e(a) ■- [ - /, . dg . (24) 

J a — e{q) + is 

We have 

Im 0(a) = —TT<P{a) (25) 

with <!> defined in (9). 

We rewrite the Hamiltonian as H — Hq + V with 

Ho ■- uj{p) := e(p) + \^e{p), H := AH - \^e{p) . (26) 

Our renormalization includes only the lowest order self-energy. This suffices on the 
time scales we consider. 


3.2 The Expansion and the Stopping Rules 

Iterating the Duhamel formula 

^ _i f ds (27) 

Jo 

gives for any fixed integer N > 1 

N-l 

in ■■= ^ V’n(i) + ^Nit) , (28) 

n=0 

with ^ 

^„(t) :=(-*)" / [dsjjyi (29) 

Jo 

being the fully expanded terms and 

(30) 

Jo 

is the non-fully expanded or error term. We used the shorthand notation 

Since each potential V in (29), (30) is a summation itself, V = —X^9{p) + A^ V, 
both of these terms in (29) and (30) are actually big summations over so-called 
elementary wavefunctions, which are characterized by their collision history, i.e. by 
a sequence of obstacles labelled by 7 e Z'^ and a label corresponding to an insertion 
oi-X^eip). 



Towards the quantum Brownian motion 


11 


Because this expansion is generated by iteration of (27), the sequences defining 
collision histories can be obtained recursively. This allows us to refine the Duhamel 
expansion by using stopping rules that depend on the type of collision history. We 
call a sequence nonrepetitive if the only repetitions in potential labels 7 occur in 
gates (immediate recollisions). The iteration of (27) is stopped when adding a new 
entry to the sequence makes it violate this condition. This can happen because of a 
recollision, a nested recollision, or a triple collision. The precise definition of these 
recollision types is given in [11]. The iteration is also stopped when the last entry 
in the sequence causes the total number of gates and i9’s to reach 2. If the sequence 
stays nonrepetitive and the total number of gates and i9’s stays below 2 , the iteration 
is stopped when the number of non-gate potential labels reaches 

K = (31) 

Note that K is much bigger than the expected typical number of collisions, 

We denote the sum of the truncated elementary non-repetitive wave functions 
with at most one power from the non-skeleton indices or ti’s and with K skeleton 
indices by V’if The superscript (< 1) refers to the number of gates and i9’s, 
each of which gives a factor A^. By this splitting, we arrive at the following modified 
Duhamel formula, in which all non-error terms are nonrepetitive. 


Proposition 1. [Duhamel formula] For any K > 1 we have 

K-l 

I —itH I \ '' / (<l),nr 

V't = e V'o = 2^ 

fc =0 


(32) 


-b f dse"''*"")" 


Ip 


(<l),nr 

*s,K 


k=0 




(<l),rec 

*s,k 


+ 


(<l),nest 

*s,k 


+ '(p 


*s,k j 


The terms under the integral correspond to the various stopping criteria indicated 
above. For the precise definition of the corresponding wave functions, see [11]. 

The main contribution comes from the non-repetitive sequences with k < K, i.e. 
from the first term in (32). The estimate of the terms in the second line (32) first 
uses the unitarity of the full evolution 


ds e 




< t ■ sup I 

S<t 


(33) 


For [[ = rec, nest, tri we will use the fact the Feynman graphs arising in the expec¬ 
tation Ejli/jf'lp contain an additional oscillatory factor, which renders them smaller 
than the corresponding non-repetitive term. It turns out that the oscillation effect 
from one single recollision, nest or triple collision is already sufficient to overcome the 
additional factor t arising from the crude bound (33). This fact relies on estimates 
on singular integrals concentrating on the energy level sets X'e. It is a well-known 
fact from harmonic analysis, that such singular integrals can more effectively be 
estimated for convex level sets. This is why the non-convexity of the energy shells 
is a major technical complication for the discrete model in comparison with the 
continuous case, where the level sets are spheres. 

Non-skeleton labels also give rise to a smallness effect due to a cancellation 
between gates and i9’s, however, one such cancellation would not be sufficient to beat 
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the t-factor. This is why at least two such cancellations are necessary in 
Finally, the term is small because it has unusually many collisions, thanks 

to the additional factor A“'* in the definition of K. 

In this exposition we focus only on the non-repetitive terms, , because 

estimating them involves the main new ideas. The error terms are estimated by 
laborious technical modifications of these ideas. 


3.3 The norm of the non-repetitive wavefunction 

We first estimate the norm of the fully expanded wave function with no gates or 
This is the core of our analysis. 


Feynman Graphs 


The wavefunction 

7 "'o 

where the summation is over all sequences for which the potential labels 7 i are all 
different. Therefore every term in 

7 , 7 ' 

has 2k potential terms, and their expectation. 


is zero, using EViy = 0, unless the potentials are paired. Since there is no repetition 
within 7 and 7', all these pairings occur between 7 and 7', therefore every pairing 
corresponds to a permutation on {1, 2,..., k}. The set of such permutations is de¬ 
noted by Vk and they can be considered as a map between the indices of the 7 and 
7 ' labels. 

We recall the following identity from Lemma 3.1 of [10] 


[ds; 


ifc+i 


fc+i 

n* 


-isjbj(pj) _ 


e!! 

27r 


da e 


k + l 

n 


i=i i=i 

for any r; > 0. We will choose 77 := Therefore, we have 

\ 2 k^ 2 tr] 

I (0) ,nr II 2 _ ^ ^ 


a — uj{pj) + irj 


Ell 


( 2 . 


E E /dpdp S{pk+i — Pk+i) 


(34) 


71,.-,7fc 

7i^7,- 


k _ _ 

XE n y-tj (Pj+i - PjWij iPcT(j)+i - Pa(j))M{k, p, p, p)ipo{pi)i>oipi) 
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with p = {pi,P2, ■ ■ .,Pk+i), f dp := f^^a)k+i dpidp 2 • • • dpfc+i, similarly for p and 
dp, and 


M^(k,p,p) ■■= 


p p / fc+l 

j J dadp . 


a — u){pj) —ipfd — oj{pj) + ip 


We compute the expectation: 


(35) 


k _ k 

ri.-'.Tfc j=i 

(36) 

In the continuous model this formula also contains a product of B-terms, where B 
was the single site potential function in (5). These factors are included into the defi¬ 
nition of Mrj ■ Most importantly, they provide the necessary decay in the momentum 
variables in the case of non-compact momentum space. A similar idea was used in 
[ 10 ]. 

Due to the restriction 7 ; 7 ^ 7 ^ , (36) is not a simple product of delta functions 
in the momenta. We have to use a connected graph expansion that is well known 
in the polymer expansions of field theory (see, e.g. [23]). We do not give the details 
here, we only note that the result is a weighted sum over partitions of the index set 
{1,..., fc}. Each term in the sum is a product of delta functions labelled by the lumps 
of the partition and each delta funtion imposes the Kirchoff Law for the incoming 
and outgoing momenta of the lump and its cr-image. The trivial partition, where 
each lump has a single element, carries the main contribution. Estimating the terms 
with nontrivial partitions can be reduced to estimates for the trivial partition [ 11 ]. 
We therefore discuss only the contribution from the trivial partition to 
given by 


V„{k,a) 


)2k^2tr, 


dpdpM^(fc, p, p)(5(pfc+i - pfe+i)^/)o(pi)^(>o(pi) 


k 

X (pi+1 - Pi - {p,T(i)+i - P<T(i))) (37) 

i=l 

This complicated formula can be encoded by a Feynman graph and Vri{k, a) is called 
the value or amplitude of the graph. The Feynman graph for the trivial partition 
corresponds to the usual Feynman graphs for the Gaussian case discussed in [10] 
and we briefly describe their construction. A Feynman graph consists of two directed 
horizontal lines (upper and lower) with k collision vertices on each that represent 
the collision histories of ip and ip, respectively. These two lines are joined at the 
two ends. This corresponds to evaluating the L^-norm on one end and inserting 
the initial wavefunction ipo on the other end. Each horizontal segment carries a 
momentum, pi,p 2 , ■ ■ -Pk+i and pi,p 2 ,. ■ -Pk+i and a corresponding (renormalized) 
propagator, (a — u{pj) — ip)~^ and {fi — u>{pj) -I- ip)~^■ Here a and /I are the dual 
variables to the time on each line and they will be integrated out, see (35). Finally, 
the collision vertices are paired. Each pairing line joins an upper and a lower vertex 
and thus can be encoded with a permutation u £ "Pa,. It is useful to think of the 
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momenta as flowing through the lines of the graph. The delta function associated 
to each pairing line in the value of the graph (37) then expresses the Kirchhoff Law 
for the flow of momenta adjacent to the two vertices. 

A typical graph with trivial partition is shown on Fig. 1. For the special case of 



Fig. 1. Typical Feynman graph with no lumps 


the identity permutation a = id we obtain the so-called ladder graph (Fig. 2). The 
following proposition shows that the ladder gives the main contribution. 






f-f-f-f- f- 


Fig. 2. Ladder graph 


The Main Contribution is the Ladder 

Proposition 2 (L^-estimate). Let r]~^ ~ t, t — 0(A^''"'‘) and k < K := 

For sufficiently small X, k and 5 there exists a positive number ci{k, 5) such that 

= K,(fc, id) -f 05 (38) 

The threshold values for k, S and the explicit form of ci(k, S) are found in [11]. 

Sketch of the proof. As mentioned above, we discuss only how to estimate the 
contributions from the trivial partition, but for an arbitrary permutation a. 

Given a permutation ct £ "Pfe, we dehne a (fe -|- 1) x (fc -|- 1) matrix M = M{a) as 
follows 

r 1 if 5-{j - 1) < i < 5-{j) 

Mij{a) := -1 if d(j) < i < 5{j - 1) (39) 

[ 0 otherwise 

where, by dehnition, a is the extension of cr to a permutation of {0,1,. .., fc -1- 1} by 
o'(O) := 0 and d{k -|- 1) ~ fc -|- 1. It is easy to check that 
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X 2 k 2 tv r s 

Vnik, a) ■- / dpdp p, p) Pi - ^ Mijpjj , (40) 

i=l 

in other words, the matrix M encodes the dependence of the p-momenta on the 
p-momenta. This rule is transparent in the graphical representation of the Feynman 
graph: the momentum pj appears in those pi’s which fall into its ’’domain of depen¬ 
dence”, i.e. the section between the image of the two endpoints of pj, and the sign 
depends on the ordering of these images (see Fig. 3) 



These momenta equal + P j+ ••• These momenta equal -p j+ ... 

Fig. 3. Domain of momenta dependencies 


The matrix M{a) has several properties that follow easily from this structure: 

Lemma 1. For any permutation a GVk the matrix M{a) is 
(i) invertible; 

(a) totally unimodular, i.e. any subdeterminant is 0 or ±1. 

The following definition is crucial. It establishes the necessary concepts to mea¬ 
sure the complexity of a permutation. 

Definition 1 (Valley, peak and slope). Given a permutation a £ Vk let a be 
its extension. A point {j,<j{j)), j £ Ik := {1, 2, ...,fc}, on the graph of a is called 
peak if a{j — 1) > a{j) < a{j + 1), it is ealled valley if a{j — 1) < a{j) > 
a{j + 1), otherwise it is called slope. Additionally, the point (fc -|- 1, fc -|- 1) is also 
called valley. The set I = {1, 2,..., fc -I- 1} is partitioned into three disjoint subsets, 
I = 7„ U /p U Is, sueh that i £ Iv,Ip or Is depending on whether {a~^{i),i) is 
a valley, peak or slope, respectively. Finally, an index i £ Iv G Is is called ladder 
index if \a~^ (i) — a~^ {i — 1)\ — 1. The set of ladder indices is denoted by It (Z I and 
their cardinality is denoted by £ = £(a) ~ \Ijg\. The number of non-ladder indices, 
d{a) := k 1 — £(a) is called the degree of the permutation a. 


Remarks: (i) The terminology of peak, valley, slope, ladder comes from the graph of 
the permutation a viewed as a function on {0,1,. .., fc -I- 1} in a coordinate system 
where the vertical axis is oriented downward. 

(ii) For a — id we have Ip = Is = {1, 2, ...,fc}, = {fc -|- 1} and It = 
{1, 2,..., fc -I- 1}. In particular, d{i(l) = 0 and d{o) > 0 for any other permutation 
a yf id. 
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The following theorem shows that the degree of the permutation d{a) measures 
the size of a). This is the key theorem in our method and we will sketch its 
proof separately in Section 3.4. 

Theorem 3. Let rj~^ := t, t = with a sufficiently small k. Let a £ Vk 

and assume that k < K = X~^{\^t). For suffieiently small k and 5 there exists 
C 2 (k,S) > 0 such that 

fr) I < A<1. (41) 

This theorem is complemented by the following lemma: 

Lemma 2. Let k = 0(A“'‘“^), d > 0 integer and let y > k + S. Then 

(42) 

cr^Vk 
<i(cr)> d 

for all sufficiently small X. 

The proof follows from the combinatorial estimate on the number of permutations 
with a given degree: 

#{a £ Vk ■■ dffi) ^d}< {Ckf . 

From Theorem 3 and Lemma 2 we immediately obtain an estimate on the contribu¬ 
tion of the trivial lumps to if s^nd 5 are sufficiently small: 

^ (43) 

cT^id 

with some appropriate C 3 {it,S) > 0. 


3.4 Sketch of the proof of the main technical theorem 

In this section we explain the proof of Theorem 3. We set 


E^{M) ■- A^ 


fc+i 


dad/l /dpH 


fc+i 

n 


|a - uj{pi) - tri\ 1 i \p _ Mjtpi) + ffil 

(44) 

For the continuous model, the definition includes the B factors to ensure the integra- 
bility for the large momentum regime. It is easy to check that Vrj{k,a) is estimated 
by Erf{M{(j)) modulo constant factors and negligible additive terms coming from 
the regime where a or /3 is big. 

The denominators in this multiple integral are almost singular in certain regimes 
of the high dimensional space of all momenta. The main contribution comes from 
the overlap of these singularities. The overlap structure is encoded in the matrix 
M, hence in the permutation a, in a very complicated entangled way. Each variable 
Pj may appear in many denominators in (44), so successive integration seems very 
difficult. We could not find the exact order (as a power of A) of this multiple integral 
but we conjecture that true order is essentially Our goal in Theorem 3 is to 
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give a weaker bound of order i.e. that is still a A-power linear in the degree, 

but the coefficient considerably smaller than 2. 

Notice that the a-denominators in (44) correspond to the columns of M and the 
/?-denominators corresponds to the rows. For this presentation we will use j to label 
row indices and i to label column indices. We recall the sets Iv,Ip, h from Definition 
1 and we will view these sets as subsets of the row indices of M. 

First we notice that if j € {le \ Iv), i.e. j is a non-valley ladder row, then there 
exists a column index i = c(j) such that the momentum pi appears only in the 
j-th /3-denominator. In other words, the i-th column of M has a single nonzero 
element (that is actually ±1) and it is in the j-th row. Therefore the dpi integral 
can be performed independently of the rest of the integrand by using the following 
elementary but quite involved bound for small k: 


sup 

w,o',f3 


dp 


T<i |a - w(pi) - ipl 1/3 - w(±pi + w) + iv\ 


< l-bO(A^/^) . 


(45) 


Note that the constant of the main term is exactly 1. This fact is important, since 
in graphs with low degree this estimate has to be raised to a power \Ie\ Iv\ that 
may be comparable with k. Clearly for k < K ^ and k -|- 5 < 1/4 we have 

-I- ©(A^^"^)^ < const , (46) 

but had 1 been replaced with a bigger constant in (45), we would obtain an expo¬ 
nentially big factor (const) ^ that would not be affordable. The precise constant 1 in 
the estimate (45) is related to the appropriate choice of the renormalization 9(p) in 
oj{p). 


After the non-valley ladder rows are integrated out, and the corresponding rows 
and columns are removed from the matrix M, we obtain a smaller matrix 
describing the remaining denominators. In we keep the original labelling of 

the rows from M. 

Now we estimate some of the /3-denominators in (44) by norm, i.e. by 
This is a major overestimate, but these denominators are chosen in such a way 
that the entangled structure imposed by M becomes much simpler and many other 
denominators can be integrated out by L^-bounds that are only logarithmic in A. 

We start with estimating all /3-denominators in rows j £ Ip by the trivial L°°- 
norm. The corresponding rows are removed from in this way we obtain a 

matrix . Let 

7* :=7\ (7pU(7,\7„)) 

be the remaining row indices after removing the peaks and the non-valley ladders. 

Then we inspect the remaining rows j £ 7* of in increasing order. The 

key observation is that for each j £ 7* there exists a column index, i = c{j), such 
that the variable pi appears only in the j-th /3-denominator, provided that all /3- 
denominators with j' < j have already been integrated out. In view of the structure 
of it means that for any j £ 7* there exists a column i = c{j) such that the 

only nonzero element among {M^p : j' > j} is M^p. This fact follows from the 
structure of M{a) and from the fact that all rows with j £ Ip have been removed. 

This property allows us to remove each remaining //-denominator, one by one, 
by estimating integrals of the type 
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1 

\a-ZJ{pi) - iri\ 


\P 


1 _ 

Lj{±pi + w) + irj\ ~ |w| 


(47) 


where w is a linear combination of momenta other than pi. The absolute value |w| 
is interpreted as the distance of w from the nearest critical point of the dispersion 
relation e{p). The variable pi at this stage of the procedure appears only in these 
two denominators. 

The exponent r can be chosen zero (with logarithmic corrections) for the con¬ 
tinuous model and this fact has already been used in [10]. For the discrete model 
we can prove (47) with r = 3/4 -|- 2 k and we know that the exponent cannot be 
better than 1/2. The reason for the weaker estimate is the lack of convexity of the 
level set Eg. Replacing a)(p) with e{p) for a moment, the inequality (47) with t = 0 
essentially states that the level set {a = e(p)} and its shifted version {(3 = e(p + w)} 
intersect each other transversally, unless w is close to zero. Indeed, the transversal 
intersection guarantees that the volume of the p values, where both denominators 
are of order rj, is of order rj^. Then a standard argument with dyadic decomposition 
gives the result with a logarithmic factor. For translates of spheres the transversal 
intersection property holds, unless w ~ 0. However, in certain points of the level 
sets Ss of the discrete dispersion relation the curvature vanishes, in fact Fe even 
contains straight lines for 2 < e < 4. The transversal intersection fails in certain 
regions and results in a weaker bound. 

Neglecting the point singularity |w|~^ in (47) for a moment (see Section 
3.5 later), we easily see that with this algorithm one can bound En{M{a)) by 
\' 2 (k-q)^-P^-T(k-p-q)^ modulo logarithmic factors, where p = |/p| is the number 
of peak indices and q ~ \Ie\ Iv\ is the number of non-valley ladder indices. From 
the definitions it follows that the sets Ip and Ie\Iv are disjoint and |/„| =p-\-l. 
Thus we have 2p-|-l-|-g<fc-|-l. Therefore 


y2(k-q) ^-p{l-T)-T{k-q) ^ ^yA_^T + l-^{k-q)/2 ^ ^ ^ [d(o-) -1] 

since q < €. If r < 1, then with a sufficiently small k we see that is a positive 

power of A. Thus we obtain a bound where the exponent of A is linear in d{a). 
With a more careful estimate one can remove the additional —1 in the exponent. In 
particular, for the continuous case with t — 0 this argument works up to k < 2 . 


We end this section with a remark. Apparently the bound k <2 (or, equivalently, 
t <C A“"^) shows up in two different contexts in this argument. To avoid misunder¬ 
standings, we explain briefly that neither of these two appearences is the genuine 
signature of the expected threshold k = 2 for our expansion method to work. The 
true reason is the one mentioned in the introduction: even the best possible bound, 

A2d(.)^ 

on the graph with permutation a, cannot beat the k\ combinatorics of the 
graphs beyond k — 2. 

In the argument above, on one hand, k < 2 is related to the error term in 
the ladder calculation (45). This error term can be improved to A^| log A| and it is 
apparently due to the fact that the renomalization term 6{p) was solved only up to 
lowest order. An improvement may be possible by including more than the lowest 
order of the self-energy. 

The second apperance of k < 2, or t ^ A“^, at least for the continuous model, 
is in (48) and it is due to the fact that certain /I-denominators are overestimated 
by L°°. This is again a weakness of our method; we did overestimates in order to 
simplify the integrand. 
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3.5 Point singularities 


The argument in the previous section has neglected the point singularity arising 
from (47). While a point singularity is integrable in d > 3 dimensions, it may 
happen that exactly the same linear combinations of the independent variables keep 
on accumulating by the repeated use of the bound (47). In that case at some point 
a high negative power of |w| needs to be integrated. While it is possible to improve 
the estimate (47) by changing the denominator on the right hand side to |w| + rj, 
this would still yield further negative jy-powers. 

It is easy to see that this phenomenon does occur. Primarily this would have 
occurred if we had not treated the ladders separately: if pi’s are ladder variables, 
then the corresponding w momenta in (47) are indeed the same. Although we have 
removed the ladders beforehand, the same phenomenon occurs in case of a graph 
which contains ladder only as a minor but not as a subgraph. Our separate ladder 
integration procedure (45) can be viewed as a very simple renormalization of the 
ladder subgraphs. The correct procedure should renormalize all ladder minors as 
well. 

To cope with this difficulty, we have to follow more precisely the point singular¬ 
ities. To this end, we define the following generalization of En{M). For any index 
set I' G I = {1, 2,... ,k + 1}, any \I'\ x (A: -|- 1) matrix M, any v integer and any 
n X (fc -I- 1) matrix £ we define 


Er,il',M,£) ~ 



dp 


( TT 1 

1 

TT 

1 11 \a-cJipi) 
\ iei' 

- 1/? - w( MijPj) + iri\ y 

IT=1 




(49) 


We follow the same procedure as described in Section 3.4, but we also keep track of 
the evolution of the point singularity matrix £. At the beginning I' = I, n = 0 and 
£ is not present. After the first non-ladder type integration, a point singularity will 
appear from (47). Some of the point singularities may get integrated out later as one 
of their variables become integration variable. Therefore we will need the following 
generalization of (47): 


Lemma 3. There exists a constant C such that for any index set A 

\cc\Tm<idj |a - w(j5) - \I3 - u{r + p) E ip\ ^ |ra-|-J3| 

< I logr;|® ^ 

a^A a '' III 

aVa 

For the eontinuous model t' = 0 while for the diserete model t' = j + 2fv. 

Using this lemma, we can keep record of the evolution of the point singularity matrix 
£ at an intermediate step of our integration algorithm. These matrices change by 
simple operations reminiscent to the Gaussian elimination. 


( J. J. Ir.. — r„/1 ) Irl 


(50) 
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Three complications occur along this procedure, we briefly describe how we 
resolve them: 


(1) The inequality (50) does not allow higher order point singularities. Although 
it is possible to generalize it to include such singularities as well, we followed a 
technically simpler path. In addition to the indices j € Ip, we select further /3- 
denominators that we estimate by the trivial bound. These additional indices 
are chosen in such a way, that (i) the number of remaining rows be at least id(cr); 
(ii) the point singularity matrix be of full rank at every step of the algorithm. This 
second criterion guarantees that no higher order point singularities occur. Since 
every £ can be derived from M by a procedure that is close to Gaussian elimination 
and M is invertible (Lemma 1), the full-rank property is relatively easy to guarantee. 

(2) The full-rank property actually needs to be guaranteed in a quantitative way, 
at least the entries of £ needs to be controlled. These entries appear in the point 
singularity denominators of (50) and their inverses would appear in the estimate. 
The key observation is that each entry of every matrix £ along the procedure is 
always 0, 1 or —1. It is actually easier to prove a stronger statement, namely that 
every f is a totally unimodular matrix. The proof follows from the fact that every £ 
can be derived from M by elementary Gaussian elimination steps plus zeroing out 
certain rows and columns. Such steps preserve total unimodularity and M is totally 
unimodular by Lemma 1. 

(3) After all /^-denominators are eliminated, we are left with an integral of the 
form 


Ep{J,th,£) 


da 


ndp)n 


n 


|a - u{pi) - i??| -li i £riPi 


(51) 


for some index set J and some point singularity matrix obtained along the integra¬ 
tion procedure. Without the point singularities, this integral could be estimated by 
the I J|-th power of | logr;|. Since £ is totally unimodular, a similar estimate can be 
obtained for (51) as well. 


4 Computation of the main term and its convergence to 
a Brownian motion 

Our goal is to compute the Wigner distribution {X, v) with t = and 

s — From Proposition 2, and similar bounds on the repetitive terms in 

(32), we can restrict our attention to the ladder graph. The following lemma is a 
more precise version of the ladder integration Lemma (45) and it is crucial to this 
computation. We present it for the more complicated discrete case. The proof is a 
tedious calculation in [11]. 

Lemma 4. Suppose f{p) is a function on T'^. Recall 0 < k < 1/16 and define 
7 := (a -I- /3)/2. Let p satisfy < 7 < Then for jrj < ipg have, 

f _ AV(p) _ 

J {a — uj{p —r) — ip){l3 — uj{p +r) + ip) 


(52) 
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= —2niX^ 


(a 


_ f(p) Sjejp) - 7 ) _ 

P) + 2(Ve)(p) • r — 2 i[A^Im 0 ( 7 ) + rj] 


dp + 0(A^/"-®"|logA|) . 


Since the Boltzmann collision kernel is uniform on the energy shell, the calcu¬ 
lation of {X, v) is more straightforward for the discrete case. We present the 

sketch of this calculation, the continuous model requires a little more effort at this 
stage. 

Let e = ^ 2 +k /2 space scale. After rescaling the Wigner function at time t, 

we compute W{e^,v) tested against a smooth, decaying function 0{^,v). In partic¬ 
ular ^ is of order 1. After the application of Lemma 4 (with v = Vk+i) and change 
of variables a ~ (a + /3)/2 and b := A“^(a — P), we obtain 

{0,mV) := f dvd^O{p,v)-EW{e^,v) = Yl [ f / dv 

k<K'^ ^ 


fc +1 


x(!?(C,Ufc+i)Wo(ee,ui)]^ 


E 


dad 6 

e 


(2 


n 


1 - oj{vi + ^) - ir]^ (^P - oj{vi - ^) + 

—2'Ki 5{e{vj) — a)dvj 


b + X~^£Ve(vj) ■ ^ — 2il(a) 


Wo(£^,vi)0(^,vk+i) 


where we defined T('y) := Im0('y) for brevity. We used rj = A^^'^ to estimate the 
error terms. The main term (left hand side above) however, is independent of rj, so 
we can choose rj = for the rest of the calculation and we note that Lemma 4 

holds for this smaller rj as well. This is the reason why the factor is negligible. 

We expand the fraction up to second order in e, we get 


—i 

b + \-^£Xe{vj) ■ ^ - 2il{a) 


—i 

b — 2iT{a) 


A-\Ve(uj)-g A-V[Ve(uj)-C]^ 
6 - 2iT{a) (6 - 2iT{a)Y 


By symmetry of the measure 27rb(e(v) — a)dv under the sign flip, v —> —v and using 
(Ve)(u) = —Ve(—w), we see that the first order term vanishes after the integration. 
We also define the matrix 


D{a) 


1 

21 ( 0 ) 


dfla{v) 


Ve(v) 

27T 


Ve(v) 

271 


After integrating out all momenta and changing the b variable we obtain 


(O, EW) « 

k<K 



2ir(a)da 


((!?(C,-))a(VLo(eC,-))a 



^2i\^thT{a) 


X 



X 


(27r)\^A-^C.^M0 

21{a) 


1 


fc + 1 


(&- + 


We sum up the geometric series and perform a residue calculation to evaluate the 
d& integral. We obtain that the main contribution comes from k ~ 2\^tT{e), so the 
truncation k < K can be neglected and the result of the dfe integration is 
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(27r)^£^f(^,£)(e)0 ^ 

To obtain a nontrivial limit, ~ 1 is necessary. Noting that t — with 

T = 0(1), we see that indeed the space must be scaled by e = Finally we 

obtain 

{0,E1T)« j Aij ^Ig^(O(^,.))a(Wh(£C,-))aexp(-(27r)^r(^,O(a)0) 

Since exp[—(27r)^r(^,Z)(a)^)] is the fundamental solution to the heat equation (15), 
from the definition of (•), and after inverse Fourier transform we obtain (17). This 
completes the sketch of the calculation of the main term. 
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